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Correlation between Interstellar Polarization and Dust Temperature: 
Alignment of Grains by Radiative Torques is Ubiquitous? 
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Abstract 

We investigate the efficiency of interstellar polarization pa/^Aj where is the fractional linear polariza- 
tion and Ax is extinction, in 14 lines of sight as a function of wavelength A. We have used the data of lines 
of sight to the Pleiades cluster obtained with the low-dispersion spectropolarimeter HBS as well as those 
in literature. It is found that the polarization efficiency p\/Ax is proportional to cxp(--/3/A) in wavelength 
A w 0.4 ^ 0.8/U,m, where /? is a parameter which varies from 0.5 to 1.2 ^m. We find that f3 is negatively 
correlated with the dust temperature deduced from infrared data by Schlegel et al., suggesting that the 
polarization efficiency is higher in short wavelength for higher temperature. According to the alignment 
theory by radiative torques (RATs), if the radiation is stronger, RATs will make small grains align better, 
and the polarization efficiency will increase in short wavelength. Our finding of the correlation between /? 
and the temperature is consistent with what is expected with the alignment mechanism by RATs. 
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1. Introduction 

Observed linear polarization in the light from distant 
stars, i.e., often called as "interstellar polarization", is in- 
terpreted phenomenon of dichroic extinction, and 
shows that interstellar grains are optically anisotropic 
and aligned, although the mechanism of the alignment 
has been still on debate (Lazarian 2007, for a recent re- 
view). The alignment of grains had been explained with 
the paramagnetic relaxation of thermally spinning grains 
that obtain angular momentum by collisions with gas par- 
ticles (Davis & Greenstein 1951, hereafter DG). However, 
the DG mechanism is not efficient, and it cannot explain 
the interstellar polarization quantitatively. For more effi- 
cient alignment, Purcell (1979) assumed a spin-up of grain 
by the ejection of molecular hydrogen from grain surface 
(the "pinwhccl mechanism"), though there still remain 
problems in quantitative explanations (Lazarian & Drainc 
1999; Roberge & Lazarian 1999). 

Dolginov & Mitrofanov (1976) first pointed out that ir- 
regularly shaped grains that have " helicity" can spin up by 
radiative torques (hereafter RATs) . More recently, Draine 
& Weingartner (1996) showed that RATs are very effective 
to align grains. Since magnetic moments within rotating 
grains are induced by the Barnctt effect, the grains preccss 
around the magnetic field, and the direction of alignment 
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is usually parallel to the interstellar magnetic field (e.g. 
Draine & Weingartner 1997; Lazarian & Hoang 2007), 
i.e., the same direction as that by the DG mechanism. The 
alignment by RATs can be more efficient if grains have su- 
perparamagnetic inclusions (Lazarian & Hoang 2008) or 
if the pinwheel mechanism is working with RATs (Hoang 
& Lazarian 2009). 

The efficiency of the RATs alignment varies with 
strength and spectral energy distribution of the radiation 
field, and thus the size of aligned grains should vary ac- 
cordingly (Draine & Weingartner 1996; Cho & Lazarian 
2007). Observationally, the maximum wavelength Amax 
of polarization in dark clouds was shown to be corre- 
lated with extinction Ay in the V-band (Whittct ct al. 
2001; Andersson & Potter 2007). Andersson & Potter 
(2010) showed that grain alignment is enhanced by the 
stellar radiation in the vicinity of a young star HD 97300 
in the Chamaeleon I cloud. For stars in the Taurus dark 
cloud, Whittet et al. (2008) showed that the polariza- 
tion efficiency pk/tk in the K-band, where tk is opti- 
cal depth, decreases smoothly with Ay beyond the re- 
gion where ice mantle feature was detected. This sug- 
gests that the alignment efficiency is not directly related 
to the state of grain surface, as is expected by the "pin- 
wheel" alignment. Those observations suggest that the 
RATs alignment works in dark clouds and star forming 
regions. However, it is still not clear whether RATs align- 
ment works or not in more diffuse clouds. 

The Pleiades cluster is associated with the diffuse re- 
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flection nebula, where grain alignment may be enhanced 
by strong stellar radiation if the alignment by RATs 
works. This motivated us to observe polarization in the 
lines of sight to stars in the Pleiades cluster with the 
low-dispersion spectropolarimeter HBS (Kawabata et al. 
1999). In this Letter, using our polarimetric data and 
those available from literature, we investigate correlations 
between polarization quantities and dust temperature, be- 
cause such correlations may be expected from the RATs 
alignment theory. 

2. Observations and Data Reduction 

We observed 8 stars in the Pleiades cluster with the low- 
dispersion spectropolarimeter HBS (Kawabata et al. 1999) 
attached to the 1.88m telescope in Okayama Astrophysical 
Observatory, from October 30 to November 13 in 2008, 
and from January 13 to 19 in 2009. HBS has a super- 
achromatic half-wave plate and a quartz WoUaston prism, 
and it can measure linear polarization over a wide range 
of the optical region, from 0.4 to 0.8/ini. We used a 
slit of 0.2mm width, yielding a wavelength resolution of 
~ 0.006/im. After binning, we have obtained spectropo- 
larimetric data with resolution of 0.02^m, and also syn- 
thesized V-band data. A unit of the observing sequence 
consists of successive integrations at 0°, 22?5, 45°, and 
67?5 position angles of the half-wave plate. The exposure 
time per object was from 1800 to 17600 sec, depending on 
the brightness of object and also on weather conditions. 

The instrumental polarization was evaluated with non- 
polarized standard stars, i.e., HD 432, HD 21447, HD 
95418, HD 198149, and HD 210027 in 2008, and HD 21447, 
HD 20630, HD 95418, HD 114710 and HD 142373 in 2009 
(Kawabata et al. 1999 and references therein). The stan- 
dard deviations of the fractional Stokes parameters for 
those standard stars were 0.02-0.07% in the synthesized 
V-band. Since these values are larger than those expected 
from photon noise, i.e. 0.01-0.03% for each measurement, 
they can be interpreted as instrumental stability in po- 
larization measurement. The photon noise for program 
stars was also small, ^ 0.01%. We thus consider that 
the accuracy in our polarimetry is mainly limited by the 
instrumental stability, and estimated the error for frac- 
tional polarization py to be 0.04%, which is a mean value 
of standard deviations for nonpolarized standard stars. 

The position angle 9v of linear polarization in the V- 
band was calibrated with polarized standard stars, i.e. HD 
7927, HD 43384, and HD 204827 in both observational 
runs. Since HD 43384 shows small temporal variation 
in polarization (Hsu & Breger 1982; Matsumura et al. 
1998), we considered it as a secondary standard, though 
no significant deviations were found from the tabulated 
value of Oy = 170° in literature. We thus calibrated 9y 
with the accuracy of ~ 0.5? 

From the observation of nonpolarized standard stars 
through a Glan- Taylor prism, we estimated the instru- 
mental depolarization to be 0.95-0.99 depending on wave- 
length, and used it for calibration. The position angle 9\ 
in wavelength A = 0.4 ~ 0.8/xm was calibrated with refer- 
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Fig. 1. Wavelength dependence of (a) fractional polarization 
px and (b) position angle 6x for the Pleiades stars. Solid 
lines in (a) show results of fitting by equation (1), and dashed 
lines in (b) the values of By. Data of each object are moved 
vertically in steps of 0.25 in (a) and 10° in (b). 

ence to 0v, with the same observation. 



3. Results and Discussion 

3.1. Fractional Polarization and Position Angle 

The observed fractional polarization p\ and position an- 
gle 9x are shown in Figure 1, except for HD 23985 and HD 
24118 which show low polarization ^ 0.1%. We assume an 
empirical formula for px by Scrkowski et al. (1975): 



Px =Pmaxexp(-/i:ln^(A/Amax)), 



(1) 



where Pmax is maximum polarization, if is a parameter 
that determines width of the curve, and A^ax is the wave- 
length at p = Pmax- Those derived values of Pmax, Knux, 
and K are tabulated in Table 1. Figure la shows that px 
is well expressed with equation (1). The position angle 9x 
is almost constant, though 9x of 19 Tau varies with A in 
short wavelength 1/A^ 2.3/im^^ (Figure lb). 

Since the scattered light is often strongly polarized, 
it may affect the polarimetric results of nebulous ob- 
jects, if it is not subtracted properly, and/or if the ob- 
jects arc surrounded by optically thick cloud, e.g. young 
stellar objects such as R Mon (Matsumura et al. 1999). 
However, the brightness of nebulosity around the stars in 
the Pleiades is not intensive compared with the stellar 
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Table 1. Polarization of Stars in the Pleiades Cluster 
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* Errors of py and pmax are estimated to be 0.04%. See text for details, 
t Errors of Amax, K, and 9v are written in the 2nd line for each object. 

t Deduced from the color excess by Crawford & Perry (1976), and its error is 0.02 mag. Other values of Ay are derived from 2MASS data, 
having errors of 0.04 mag. 



light, typically ^ 20 mag/arcsec^ in the B or V-band. We 
thus expect that the effect of the nebulosity is subtracted 
with the sky background in the reduction process. In ad- 
dition, the spectral types of the stars are normal. Thus the 
observed polarization is expected to be mainly foreground 
interstellar in origin. Nevertheless, 19 Tau may have a 
component of non-interstellar origin, because the position 
angle is variable: Oy was 114° in Markkanen (1977) and in 
Breger (1986), but it is 140 ±4° in our observation (Table 
1). We thus exclude 19 Tau in the following discussion. 

3.2. Polarization Efficiency 

The ^v-dcpcndence of Amax in dark clouds suggests 
that the alignment of grains is induced by RATs (see 
Section 1). However, Andersson & Potter (2007) noticed 
that the extrapolated value of Amax at Ay = in each 
cloud is correlated with the mean value of the ratio of to- 
tal to selective extinction Ry, where Ry = Ay / Eb-v and 
Eb-v is the color excess for B — V. Since Ry characterizes 
the extinction, this correlation means that Amax depends 
not only on the alignment, but also on the size of total, 
i.e. aligned and nonaligned, grains. We thus use another 
quantity less affected by the variation of grain size. 

Compared with Amax, the polarization efficiency p\/A\ 
should be less dependent on the variation of grain size, 
because such variation will be canceled in px/Ax. We 
thus explore the observed properties of px/Ax, expecting 
to obtain information on alignment. It should be noted, 
however, that Voshchinnikov & Das (2008) and Das ct al. 
(2010) showed that px/Ax depends on grain size, shape, 
material, and other parameters. It would be possible to 
examine the properties oipx/Ax in more detail with using 
light scattering calculations (e.g. Matsumura & Seki 1991, 
1996; Matsumura & Bastien 2009), but it is beyond the 
scope of this Letter. 

To evaluate extinction Ax, we have used two methods: 

Method 1: On the assumption that the A-depcndcncc 



of Ax is determined by Ry and scaled by Ay (Cardelli 
et al. 1989), we evaluate Ax with interpolating the data 
of Ax/EB-y for Ry = 2.1 5.5 tabulated in Fitzpatrick 
(2004). The values of and Ry are calculated by 
the formulae Ay = l.tSy-K and Ry = 1.1£'v-k/-E'b-v 
(Whittet & van Breda 1980), respectively, where -Ey-K 
is the color excess for V — K. We use B and V mag- 
nitudes in the Simbad database, while for the K band, 
we transform Kg magnitude in the Two Micron All Sky 
Survey (2MASS) into the K magnitude in the system of 
Koornncef (1983) with a formula by Carpenter (2001). 
We refer to Fitzgerald (1970) and Koornneef (1983) for 
the intrinsic colors of B — V and V — K, respectively. 

Since the errors of Ry derived with Method 1 are large 
for some stars (Table 1), we also use another method as 
below (Method 2). The error of Kg for 27 Tau was par- 
ticularly large, ^ 0.3 mag., we could not obtain reliable 
results, and excluded 27 Tau in the following discussion. 

Method 2: We assume that the extinction properties 
are not variable within the Pleiades cluster, and apply the 
extinction curve for HD 23512 to other stars, scaling it by 
the value of Ay deduced with Method 1. The extinction 
curve for HD 23152 is reduced by Fitzpatrick & Massa 
(2007), and most reliable among the Pleiades stars. 

We explore not only the Pleiades stars, but also the 
stars for which polarization and extinction data are avail- 
able from literature. We have used the polarimetric data 
of various stars by Weitenbeck (1999, 2004). The data of 
HD 29647 (Whittet et al. 2001) and HD 38087 (Serkowski 
et al. 1975) are used in addition to those from Weitenbeck 
(1999). Extinction data for those stars are cited from 
Fitzpatrick & Massa (2007). Also used are the data for 
high latitude clouds MBM 30 and MBM 20 (LDN 1642) 
by Seki & Matsumura (1996). 

Figure 2 shows the A-dependence ofpx/Ax, which is de- 
rived with Method 1. The values of log(pAMA) decrease 
linearly with inverse wavelength 1/A, though slight devi- 
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Table 2. Polarization Properties deduced from Literature. 
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(2.74±.54, 1.15±.12) for SAO 149760. 
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Fig. 2. Wavelength dependence of polarization efficiency 
Px/Ax deduced with Method 1 (See text). Asterisks show 
the results of HD 38087 (Serkowski et al. 1975) and filled tri- 
angles HD 29647 (Whittet et al. 2001), while other symbols 
are for the Pleiades stars as in Figure 1. Long dashed lines 
axe the results by Weitenbeck (1999). The data are normal- 
ized by a (see equation (2)), and moved vertically in steps of 
factor 1.5. 

ations from the linear relation are found. We thus make 
linear fitting in A = 0.4 ^ 0.8/im with the equation: 

In(pAMA) = lna-/3(l/A-l/0.55Mm), (2) 

where A is in /im, and a and /3 are parameters and tabu- 
lated in Tables 1 and 2. 



3.3. Polarization Efficiency and Dust Temperature 

To discuss the correlations between the polarization 
properties and dust temperature, we use the temperature 
Tdust by Schlegel et al. (1998). They deduced Tdust from 
COBE and IRAS data, on the assumption of A~^ emissiv- 
ity of grains in the infrared. Their data is homogeneous 
all over the sky, and thus suitable for our study that con- 
tains not only the lines of sight to Pleiades stars, but also 
those to other objects. 

Between Amax and Tdust, we find a weak correlation in 
Figure 3a, and the correlation coefficient r is —0.30. The 
correlations between (3 and Tdust are much better, i.e., 
r = -0.54 (Method 1, Figure 3b) and r = -0.57 (Method 
2, Figure 3c). It is remarkable that the relative positions of 
HD 210121 and HD 38087 are different between in Figure 
3a and in Figure 3bc. This is caused by the different values 
of iT!v, i-e., R\ of HD 210121 is smaller (=2.0, Fitzpatrick 
& Massa 2007), and that of HD 38087 is larger (= 5.8, 
Fitzpatrick & Massa 2007) than other objects (i?v ^ 3). 

We finally discuss the above mentioned correlations on 
the basis of the RATs alignment theory. Using their equa- 
tion (5) in Cho & Lazarian (2007) and equation (67) 
in Draine & Weingartner (1996), with typical values for 
physical quantities in interstellar space tabulated in Table 
2 of Lazarian & Hoang (2007), we can express the smallest 
size aiowcr of aligned grains as 

(aiowcr/lA*m) = 0.089 x (rdust/18K)-2, (3) 

where we assume that the efficiency Qr for RATs is 0.1, 
and that the emissivity of grains is oc A^^. 

On the size distribution of aligned grains, Mathis (1986) 
showed that the observed polarization p\ can be repro- 
duced if the fraction (1 — cxp(— (a/a')'^)) of grains with 
radius a arc aligned, where a' is a parameter for typical 
size of smallest aligned grains. Mathis (1986) then ob- 
tained the equation: 

{a'/\^lm) = 0.327 x (A„,ax/lMm)2 ". (4) 

If we assume a' = aiower (or a' = 2aiower), we can relate 
Amax and Tdust with equations (3) and (4), and draw the 
short dashed (or long dashed) line in Figure 3a. 
For the relation between j3 and Amax, we obtain 

(/3/lAim) = 1.79 X (Amax/l/^m)i-39, (5) 

with using equation (1) and extinction curve for Ry = 3.1 
(Fitzpatrick 2004). We then write /3 as a function of Tiust, 
and draw short and long dashed lines, for a' = aiowcr and 
a' = 2aiowor, respectively, in Figures 3b and 3c. 

Those lines in Figure 3 seem to follow the observations 
well, i.e., the correlation between Amax and Tdust, and that 
between /3 and Tdust can be explained by the RATs align- 
ment theory. Since those data contain regions of various 
temperature, i.e., the Taurus dark cloud, reflection nebu- 
lae, etc., our results suggest that the alignment by RATs 
is ubiquitous in the interstellar space. 
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Program at Okayama Astrophysical Observatory of 



Interstellar Polarization and Temperature 



5 



0.8 







29647 



(a) 



23512 ~~ - ^ 
• 38087 




16 18 20 

Tdust [K] 



22 



Fig. 3. (a) Correlation between 'Tdust And Amax- Filled 
circles are the Pleiades stars, open circles the results by 
Weitenbeck (1999, 2004), filled squares those by Whittet et 
al. (2001), and filled triangles those by Seki & Matsumura 
(1996). Double circle is the average of interstellar medium. 
Result of linear fitting for 14 stars is drawn as solid line. Short 
dashed line presents a model prediction for a' = aiowor i while 
long dashed line a' = 2ai^„^^ . See text for details, (b) Same 
as (a) but for T^ust and /3 with Method 1. (c) Same as (b) 
but with Method 2. 
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